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Abstract. A generalization of the eight vertex model by means of higher spin 

I representations of the Sklyanin algebra is investigated by the quantum inverse 

' scattering method and the algebraic Bethe Ansatz. Under the well-known 

ly^ , string hypothesis low-lying excited states are considered and scattering phase 

0^ ' shifts of two physical particles are calculated. The S matrix of two particle 

""^j^l states is shown to be proportional to the Baxter's elliptic _R matrix with a 

• different elliptic modulus from the original one. 

> 



Introduction 



^ ■ In this paper we consider a generahzation of the eight vertex model by means of 

. higher spin representations (spin i) of the Sklyanin algebra [Ski] on a space of theta 

functions [Sk2]. This model has 2i+l dimensional state space on each vertical edge 
and two dimensional state space on each horizontal edge. 

Relation of the eight vertex model to the SOS type model was established by 
Baxter [Bax3]. A similar relation also holds in our case and, using this relation, we 
can pursue the quantum inverse scattering method and the algebraic Bethe Ansatz, 
following [KS], [TFl], [TF2]. In the first part of this paper, we examine Bethe 
vectors and give coordinate expression of them in terms of Boltzmann weights of 
SOS type model. We also prove a sum rule of rapidities of quasi-particles, which was 
proved for the eight vertex model by Baxter [Baxl] who made use of a functional 
equation, an alternative to Bethe Ansatz. This rule is related to a parity of Bethe 
vectors. 

A higher spin version of SOS type model was constructed by Date, Jimbo, Miwa 
and Okado [DJMO], [JMO], [DJKMO] by fusion procedure [KRS], [Ch], [ZH], [HZ]. 
Recently Hasegawa showed [Ha2] that a representation of the Sklyanin algebra ob- 
tained by fusion procedure repeated 2i — 1 times is equivalent to the spin i rep- 
resentation on a space of theta functions [Sk2]. Hence, in principle our model is 
equivalent to the higher spin SOS model by Date and others. The use of repre- 
sentations by Sklyanin [Sk2] makes it possible to compute eigenvectors of transfer 
matrices explicitly and to apply the quantum inverse scattering method and the 

QlrrtiV»-rQi/^ RtifVio AnGoff? rlirtir'+Kr f 



2 



TAKASHI TAKEBE 



1 dimensional local quantum spaces. But in general the transfer matrix of our 
lattice model does not give local Hamiltonians directly, because the dimension of the 
auxiliary space is fixed to two. In order to write down the Hamiltonian of this spin 
chain, we must use the fused transfer matrix corresponding to a 2£ + 1 dimensional 
auxiliary space. More generally we can construct a model with arbitrary spins on 
quantum and auxiliary spaces by fusion procedure. We will study such models in 
the forthcoming paper. 

Note that Bethe vectors constructed above give eigenvectors of these models 
simultaneously under assumption of non-degeneracy, since transfer matrices with 
different auxiliary spins are mutually commuting [KRS]. 

Though momenta and Hamiltonians of spin chains are calculated from transfer 
matrices of fused models, S matrices (phase shifts) of spin waves do not depend on 
the auxiliary space. In the second part of this paper, we calculate a two particle 
S matrix of spin waves from Bethe vectors obtained above, following the recipe by 
Korepin [K], Destri and Lowenstein [DL]. The result confirms Smirnov's conjecture 
[FIJKMY] which states that this S matrix should be given by an elliptic R matrix, 
the elliptic modulus of which is different from that of the original R matrix in the 
definition of the model. 

Corresponding results were established for the totally isotropic models (the 
XXX model) and its higher spin generalization, by Faddeev, Takhtajan, Babu- 
jian, Avdeev and Dorfel [TF2], [Takh], [Bab], [AD] and for the XX Z model (the 
six vertex model) and its higher spin generalization, by Sogo, Kirillov, Reshetikhin 
[So], [KiR]. Free energy of the eight vertex model was obtained by Baxter [Baxl] 
and low- lying excited states were studied by Johnson, Krinsky and McCoy [JKM], 
but our calculation of the S matrix seems to be new even for the eight vertex model 
{£ — 1/2), though a partial result on the S matrix for this case was calculated by 
Freund and Zabrodin [FZ]. The algebraic Bethe Ansatz was shown to be appli- 
cable to the higher spin eight vertex models in [Takel] and their free energy was 
calculated in [Take2] . 

This paper is organized as follows: In Chapter I we begin with review of def- 
inition of the model and generalized algebraic Bethe Ansatz, following [Takel]. 
Then, giving a whole set of intertwining vectors explicitly, we write down coor- 
dinate expression of Bethe vectors in terms of them. A sum rule of rapidities of 
quasi-particles arc presented which helps solving the Bethe equations. The proof is 
given in Appendix B. In Chapter II we study thermodynamic limit of several Bethe 
vectors. Free energy is calculated in §2.2 (this result was announced in [Take2]) and 
low-lying excited states are examined in §2.3 under assumptions of string configua- 
tions. In particular we compute a two particle S matrix in §2.3. We summarize 
prerequisites on the Sklyanin algebra and its representations in Appendix A. 

Chapter I and Appendices are of algebraic nature, while in Chapter II we do 
not try to give mathematically rigorous arguments to every detail, since the goal of 
this chapter is to compute quantities of physical importance. In order to make this 
computation rigorous hard analysis is indispensable, which is beyond the scope of 
this paper. 
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wc change notations and normalizations. In §1.2 wc find an explicit form of a whole 
set of intertwining vectors which enables us to write down coordinate expression of 
Bathe vectors (§1.3). Only the "highest" ones of intertwining vectors ("local pseudo 
vacua" in the context of the algebraic Bethe Ansatz [TFl]) were used in [Takel] to 
construct Bethe vectors. In §1.4 we show that sum of rapidities of quasi-particles 
should satisfy an integrality condition. This comes from quasi-periodicity of theta 
functions and therefore is absent in the case of models associated to trigonometric 
and rational R matrices. 

1.1. Definition of the model. The model is parametrized by a half integer £ and 
two complex parameters: an elliptic modulus r and an anisotropy parameter rj. In 
this paper we assume that the elliptic modulus is a pure imaginary number while 
the anisotropy parameter is a rational number: 

t r 

where t > and r, r' are integers mutually coprime. Moreover we impose a condi- 
tion that r is even, r' is odd, and 2{2i + 1)?7 < 1. 

Now we define a lattice model of vertex type as in [Takel]. We consider a square 
lattice with N columns and N' rows on a torus, i.e., periodic boundary condition 
imposed. States on the n-th vertical edge belong to the spin £ representation space 
Vn ~ Qqq^ of the Sklyanin algebra (see Appendix A) while states on each horizontal 
edge are two dimensional vectors. A row-to-row transfer matrix, T(A), of the model 
is defined as the trace of a monodromy matrix, T{X), in the context of the quantum 
inverse scattering method [TFl]: 

^(^)= (c^Ia] Zlx)) ■■=^r.iX)...L.iX)L.iX), (1.1.2) 
T(A) = Trc2(r(A)) = An{X) + Dn{X), (1.1.3) 
where the L operators, L„(A), are defined by (cf. (A. 4)) 

^••<^> = few t(A)) t WiWA(S'} « 

p^(5") = 1 ® . . . ® 1 ® p^S") ® 1 ® . . . ® 1, 

elements of which act on a Hilbert space TC = (H)^=i Ki, but non-trivially only on 
the n-th component. Assignment of Boltzmann weights to vertices are determined 
by this L operator. The monodromy matrix is a 2 x 2 matrix with elements in 
Endc(H), and the transfer matrix is an operator in Endc(H). 

The partition function, Z{X), and the free energy per site, /(A), are 

Z(A) = T¥h(T(A)^'), 
-Pm = J^^ogZ{X). 
In the thermodynamic limit, A^, N' — > oo, only the greatest eigenvalue, Amax, of 



(1.1.4) 
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which was computed in [Take2]. We wiU recall this result in §2.2 with details 
omitted in [Take2]. 

Remark 1.1.1. The above defined model is a homogeneous lattice in the sense that 
it is invariant with respect to vertical and horizontal translation. We can also define 
an inhomogeneous lattice by assigning different spectral parameter, Aj, and different 
spin, £j, to each vertical edge. We have only to replace I/„(A) in (1.1.2) and (1.1.3) 

by 

3 

Ln{X - Xn) := ' An)p^(5") ® 

a=0 

P^niS") = 1® p^^iS'') ® 1 ® . . . ® 1. 

All arguments in §1.2, §1.3 remain true with suitable changes, as is shown in [Takel]. 
Such models are important for the study of certain integrable systems [NRK] . 

1.2. Intertwining vectors, gauge transformation. Intertwining vectors were 
first introduced by Baxter [Bax3], and given an interpretation as a gauge trans- 
formation in the context of the quantum inverse scattering method by Takhtajan 
and Faddeev [TFl]. Generalization to higher spin case by means of fusion proce- 
dure was studied by Date, Jimbo, Miwa, Kuniba and Okado [DJMO], [DJKMO], 
[JMO]. Here we define intertwining vectors directly in the space of theta fucntions. 
They should be identified with those defined in [DJMO], [DJKMO], [JMO] through 
Hasegawa's isomorphism [Ha2]. 

Definition 1.2.1. Let A;, k' be integers satisfying k — k' E {—2£, —2£ + 2, . . . ,2£ — 

2,2£}, and A, s = (s+,s_) be complex parameters. We call the following vectors 
0A;,fe'(A; s) = 0^^]j/(A; s){z) e Qqq'^ intertwining vectors of spin £: 



4'l'(A;«)W = aM' ri e(z+'-±^ + ^- + {k'-£ + 2j-l)ri\ 

.7 = 1 ^ ^ 



X 



k-k' 
^ 2 



•xd[z- - ^ - (fc' - £ + 2j - 1)?7 



ri ^(^ + ^^ + ^ + (/c-^ + 2j-l)r7) X 



•xe[z-'-^-\-{k-£^22-\)r) 

Here d{z) = 6'oo(-2;t) (see (A.l)), and ak,k' = e2^^^('=+'=')''(it-V2e^^(«+-s-))^, 
Following [TFl], we introduce a matrix of gauge transformation Mj.: 

M (\ N -A + 2/c?7);2it) eii{-it{s- + X + 2kr]);2it)\ 

' ^ \eQi{-it{s+-X + 2kr]y,2it) eoi{-it{s_ + X + 2kr]y,2it)J 
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where 



Wk = 2 H2A;r/- -. 



Let us define a twisted L operator by 
-t^fe,fe'(A, s) — 



■.= Mj;\X;s)L{X)Mk'iX;s). 



(1.2.3) 



(1.2.4) 



Proposition 1.2.2. Each component of L^^^i acts on the intertwining vector as 
follows: 
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k'-l 
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k' + 1 
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k'-l 
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k' + 1 



-l,fc'-l, 



'k-l,k' + l^ 



(1.2.5) 



-i,fe'-i, 



,fc'+i, 



where (j)k,k' = 0fc^i'(O)S) c-i^d W is the Boltzmann weight of SOS type [Bax3], 
[DJMO]; 
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k + 1 


k' + l 



= 2en{{k -k'- 2e)r))euiw^k+k')/2 - A), (1-2-6) 
x)^2en{X-ik-k')v)^''^''^'^''-''^/'^ 



^11 K) 

This proposition is proved in the same way as (3.7) of [Takel]. 



If we denote components of Lk^k' by 

Lk,k'{\ s) - 
(1.2.5) takes the form: 



Lk,k'{-'^, -1; A; s) Lfc,fc/(-l, +1; A; s) 
Lk,k'{+^, -1; A; s) Lk,k'{+1, +1; A; s) 



Lk,k' (s, s'; A; s)(f)k,k' = W 



k k' 
k + e k' + e' 



A (f)k+s,k'+e'- 



(1.2.7) 



(1.2.8) 



In [Takel], vectors cu^ = 0n+2^m,n+2^(m-i)(s) were called local vacua. For these 
vectors the formulae (1.2.5) reduce to 

ak,k-2e(f>k,k-2i{s) = 26'ii(A + 2ir])(j)k-i,k-2£-i{s), (1.2.9) 
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(1.3.1) 



Remark 1.2.3. Denoting the column vectors of Mk by 'i/'A;,fe±i(A; s), = 
(•i/'A;,fc-i(A; s), '0fc,fc+i(A; s)), one can rewrite (1.2.5) as follows: 

L(A - //)<^fc^i., (A; s) ipk',k'+e' (/w; s) = 

= 5Z^(a; + £ /c' + £' ^-A*) </'i+£,fc'+£'(A;s)®V'fc,fe+£(/^;«)• 
Namely (j) and intertwine the vertex weights and the SOS weights. This is where 
the name "intertwining vector" comes from. See [BaxS], [DJMO], [Hal]. Note 
that (p^i^l^i (A; s) are proportional to column vectors of under the identification 
(A.9). ' 

1.3. Generalized algebraic Bethe Ansatz. In this section we recall the con- 
struction of eigenvectors of the transfer matrix by means of the algebraic Bethe 
Ansatz, following [Takel], and give several properties of them. Hereafter we as- 
sume that M :— N£ is an integer. 

First introduce a modified monodromy matrix twisted by gauge transformation: 

T (\. ^\ - f^k,k'i^'^^) Sfc,fe/(A;s)\ 
:=M-\X;s)T{X)Mk'{X;s), 

and fundamental vectors in 7i: 

|ajv,aAr-l, • • -,01,00) := 0ajv,ajv-i <8) • • • <H) 0o2,ai <8) 0ai,ao) 

(1.3.2) 

where (f)a,b — 00,6(0; s) are intertwining vectors defined by (1.2.1). We fix a value 
of the parameter s = (s_(_,s_) and suppress it unless it is necessary. A pseudo 
vacuum, O^, is a fundamental vector characterized by gq = a, ai — Oj-i = 2£ for 
alH = 1,...,A^: 

n% = \a + 2Ne, a + 2{N - !)£, . . . , a + 2£, a). (1.3.3) 
This vector satisfies 

Aa+2Ne,ai^)n% = (2^n(A + 2erj))''Q%-\ (1.3.4) 
Da+2Ne,aWn% = (2^ii(A - 2£r])f n''+\ (1.3.5) 

Co+2Ar£,a(A)O^ = (1.3.6) 

by virtue of (1.2.9), (1.2.10), (1.2.11), respectively. 

As is shown in [Takel], the algebraic Bethe Ansatz for our case leads to the 
following: 

Proposition 1.3.1. Let v he an integer, Ai,...,Am complex numbers. Define a 
vector ^'i^(Ai, . . . , Am) by 

r-l 

*.(Ai, . . . , Am) := V e^^^'^'^'^aiXi, ■ ■ ■ , Am), 
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Then \E'jy(Ai, . . . , Am) is an eigenvector of the transfer matrix T{\) with an eigen- 
value 

t{X) = e^-^-n2^n (A + 2£r;; r))^ ^'f ~ ~ 

fj{ fiilA-AjSr) 

provided that u and {Ai, . . . , Am} satisfy the following Bethe equations; 

^n(A,- - 2£n;T)J ^n(A,- - A^ - 2rj;Ty 3 g. 

k^j y ■ ■ J 

for a// j = 1, . . . , M. 

Proof. The proof is the same as that in [TFl]. We only recall the periodicity of 
the vector with respect to a, which is the reason that we do not have to take 
an infinite sum in the definition of \E'jy. 

Recall that r] is a rational number r'/r. Therefore 2{k + r)ri = 2kri (mod 2). 
This fact and quasi periodicity of theta functions imply Mfc+^ = M^. Hence (see 
(1.3.1)) 

%+r,k'+r{X) =%,k'{X)^ (1.3.10) 

in particular, Sfc+r,fc'+r(A) = Bk,k'{>^)- Similarly one can prove (t)k+r,k'+r = (t>k,k'- 
This proves $a+r(Ai, . . . , Am) = $a(Ai, • • • , Am)- □ 

The eigenvalue (1.3.8) is written in a compact form in terms of a function Q{X) 
defined by 

M 

Q(A) = e--^ J] ^"(^ - ^i)- (1-3-11) 
The eigenvalue of the transfer matrix for a Bethe vector ^'i/(Ai, . . . , Am) is 

Ul{A) (1.3.12) 

where h{z) = (2^11(2))^. The Bethe equations (1.3.9) can be interpreted as the 
condition of cancellation of poles at \j of the right hand side of the above equation. 
This observation is due to Baxter [Baxl] and a starting point of Reshetikhin's 
analytic Bethe Ansatz [Rl] . We essentially use this observation to derive the sum 
rule in §1.4. 

Because of the commutation relation Bk.k'+i{X)Bk+i^k'{p) = Bk.k'+i{fj)Bk+i^k'{X), 
which is a consequence of relation (A. 3), ^'jy(Ai, . . . , Am) does not depend on the 
order of parameters Ai, . . . , Am- Moreover, we may restrict the parameters to the 
fundamental domain 
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Lemma 1.3.2. Suppose (z^, {Ai, . . . , Am}) is a solution of the Bethe equations. 
Then for any j, 1 ^ j ^ M , 

i) ^i,{Xi, . . . , Xj + 1,...,Xm) is proportional to ^i,{Xi, . . . , Xj, . . . , Xm), and 
{u, {Ai, . . . , Xj + 1, . . . , Am}) is a solution of the Bethe equations. 

ii) \l/;^_(_2(Ai, . . . , A j + T, . . . , Am) is proportional to \l/;y(Ai, . . . , Aj, . . . , Am), and 
{v + 2, {Ai, . . . , Aj + T, . . . , Am}) is a solution of the Bethe equations. 

Proof. The quasi-periodicity of theta functions implies that (z^, {Ai, . . . , A^ + 
1, . . . , Am}) and {v + 2, {Ai, . . . , A^ • + r, . . . , Am}) are solutions of the Bethe equa- 
tions. 

Substituting A — > A + 1 and A ^ A + r in (1.1.4), we obtain 

L„(A + 1) = ~a^L^{X)a\ 

Ln{X + T) = -e-^'^-^^'^a^Ln{X)a\ 

In the same way (see (1.2.2)): 

Mk{X+l-s)=a^Mk{X;s), 

_^ni{s++2kr]-\-r) g \ 

_g— 7ri(s_+2fe77+A) j • 

Hence from (1.3.1) follows 

Tk,k'iX+l;s) = Tk,k'iX), (1.3.14) 

/„—TTi(s++2kr] — X—T) n \ 

Tfc,fe'(A + T,s) =(-1) e I g7ri(s_+2fc^+A) ) X 

/g7ri(s++2A;'»7-A-T) q \ /i Q 1 r;\ 

x%,k' (A; s) Q ^-niis.+2k'v+x) J ■ ^^■'^■^^^ 

The (l,2)-component of (1.3.14) is Bk^k'iX + 1) = (-1)^ Bk,k'iX). Thus 
*^(Ai, . . . , A, + 1, . . . , Am) = (-l)^*^(Ai, ...,Xj,..., Am), 
which proves i). The (l,2)-component of (1.3.15) is 

Bk,k'iX + T)= const. e-''^'^>'+'''^''Bk,k'iX). 
Here const, does not depend on k, k' , but depends on A, s±. Thus 

$a(Ai, ...,Xj + T,..., Am) = const. e"^''^''''$a(Ai, . . . , A^, . . . , Am), 

and 

r-l 

*,+2(Ai, . . . , A,- + T, . . . , Am) = const. ^2.ii.+2)ar,-4.iari^^^^^^ ^ _ _ ^ 

a=0 

= const. \E'j,(Ai,A2,...,Am)- 
This proves ii). □ 



Mk{X + T;s) = a^Mk{X; s) 
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Proposition 1.3.3. Let \E'iy(Ai, . . . , Am) be as defined by (1.3.7). Then 



r-l 



*.(Ai,---,AM)-J]e2-^'^^«x 



0=0 



ao=a,ai,... \fa^.}^=lj=l \ ^ '-'J * ^.J ^ / / 
...,ajv-i,ajv=a \'- / 

Here the sum in () is taken over a set of integers aij (0 ^ i ^ M, ^ j ^ N) 
satisfying the admissibility condition 

aij - ai_ij = ±1, aij - j_i e {-21, -2^ + 2, . . . , 2^ - 2, 21}, 

and the boundary condition, 

aoj = aj, aM,j = a- N£-\- 21 j, 

aifi = a-i, ai^N = ct + i. 

Note that the Bethe equations are not assumed here. 

Proof. Operator -Bajv.ao W (1' 2)-element of the monodromy matrix 7^jv,ao W - 

LaN,aN-i ■ ■ - Lauao- Heuce by (1.2.8): 

BaN,ao{^)\(^N, OAT-I, • • • , Ol, ^o) = 

= ^ Laf^,aN-i{^N,£N-l',^) ■ ■ ■Laj,aj-i{^j,^j-l',^) ■ ■ ■ 



SN = —,eN-l,- 

...,ei,eo=+ 



(£l,£o; A) • 

AT 



£iv — — ,eiv-i,---i=l 
...,ei,eo=+ 



0a 



A'+ejv,aAr-i+ejv-i 



+ Sj Clj — 1 + — 1 



A 



''oi+ei,ao+eo' 



Applying this formula iteratively, we arrive at (1.3.16). 



□ 



1.4. Sum rule. In the previous section we defined Bethe vectors by (1.3.7). Here 
we show an integrality condition of sum of parameters A^. 

Theorem 1.4.1. Let (v, {Ai, . . . , Am}) be a solution of the Bethe equations (1.3.9). 
We assume that {Xj} satisfy the following additional conditions: For any j = 
1,...,M, 

i) Xj ^ {2{n + e)ri\ne Z}. 

a) there exists a e Z such that Xj+2ar] ^ A^ (mod Z+Zr) for any k = 1, . . . , M. 
Hi) (Technical assumption of non-degeneracy: see Appendix B.) 
Then there exist integers no, n\ which satisfy 

M 

2 ^ Aj = no + niT. (1-4.1) 

T'Vita T^T" r\r\f iG f fir»Vi n 1 r» Q 1 qti/^ /^r^T^^"t^ in ti/^ in A T^T^^in/^ i v R ATr^^"ti c» can m T^^" i r\n a i ^ 
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Baxter derived this rule in [Baxl] , directly constructing an operator on Ti which 
gives Q(A) (see (1.3.11)) as its eigenvalue. Unfortunately we have not yet found 
such an operator in our context. (Kulish and Reshetikhin [KuR] found for a rational 
R matrix case that transfer matrix "converges" to the Q operator by iteration of 
fusion procedures.) 

In addition, Baxter's result tells us that X]j=i -^j related to parities of Bethe 
vectors. These parities are associated to reversing the arrows (o"^) or to assigning 
— 1 to the down arrows {a^) of the XYZ spin chain model. They correspond to the 
unitary operators U\ and acting on the spin £ representations. (See Appendix A.) 

Lemma 1.4.2. For a = 1, 2, 3, Ua commutes with the L operator as 

U-'L{X)Ua = {a'^r^ L{X)a\ (1.4.2) 

They commute with the transfer matrix: U~^T{\)Ua = T(X). 

Proof. Adjoint action by Ua induces an automorphism Xa on the Sklyanin algebra 
(see Appendix A): 

U-'L{X)Ua^XMX)) 

= Wo{X)p\S'^) ® + Wa{X)p\S'') (g) a" 
- Wb{X)p\S') (g) - Wc{X)p\S'') (8) 

where {a,b,c) is a cyclic permutation of (1,2,3). By the anti-commutativity of 
the Pauli matrices, the right hand side of the above equation is nothing but 

((T")-iL(A)a". 

Commutativity T{X)Ua — UaT{X) is a direct consequence of (1.4.2) and (1.1.3). 

□ 

Operators on Ti. are involutive and commute with each other. In fact, by 

virtue of relations = (—1)^^ and UaUb = {—l)'^^UbUa, we have 

(urr = i-^r' = 1, 

(Recall that Nl is an integer.) Therefore an eigenvalue of Uf^ is either +1 or —1. 
Assume that \E' G is a common eigenvector of T{X) and C/a's. We assign parities 
v" and v' to ^ by 

U®N^ ^ (-1)""*, U®^^ = (-1)"'*. 

From Baxter's result [Baxl] and Theorem 1.4.1, the following conjecture seems 
to be plausible. 

Conjecture 1.4.3. Let {v, {Ai, . . . , Am}) be a solution of the Bethe equations, and 
v" and v' he parities of the Bethe vector ^'j^(Ai, . . . , Am) defined above. Then 

zy + zy' + iV£ = (mod 2), 
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2. Thermodynamic limit 

In this chapter we consider the hmit N ^ oo. Our calculation is based on 
the string hypothesis introduced by Takahashi and Suzuki [TS] which we review 
in §2.1. In §2.2 we compute the free energy of the model. In §2.3 we introduce 
four kinds of perturbation of the string configuration of the ground state. Each 
of them are parametrized by two continuous parameters which are regarded as 
rapidities of physical particles. We compute polarization of the Dirac sea of quasi- 
particles induced by these perturbation, following the recipe by Johnson, Krinsky 
and McCoy [JKM]. We also calculate eigenvalues of the S matrix of two physical 
particles, using the method developed by Korepin [K] , Destri and Lowenstein [DL] . 
The result coincides with Smirnov's conjecture [FIJKMY] . 

2.1. String hypothesis. First let us rescale the parameters so that integrals in 
later sections are taken over a segment in the real line. We denote xj = itXj. Then 
the Bethe equations (1.3.9) takes the form 

Oiiix j + 2iir)t;it) \^ ^ ^-4nir,{u+2Z^^^Xk) TT - Xk + 2ir]t;it) 

Oii{xj - 2i£r]t; it) J ^J^ Oii{xj - Xk - 2ir]t; ity^^ ^ -^^ 



while the corresponding eigenvalue is 

t{x) = (-2iVt)^e"i^^^' -4^(^+1)^'*') X 



M 



^2..,(.+2E- , ^.)^2e,,{x + 2^er}t■, ^t))^ T] ^^i^ ' ' ^ivtf) 

fj^ Oii{x-Xj;it) 

+e-'^M^+'^-^^^\2e,,{x - 2iir^t■^t))^ T] (x - + 2zryMt) 

fj^ Oii{x-Xj;it) 

where x = itX. According to Lemma 1.3.2, we may assume that | Re{xj)\ ^1/2 for 
any j. 

Now the string hypothesis [TS] is stated in the following way. For sufficiently 
large solutions of (2.1.1) cluster into groups known as ^-strings, A = 1,2, . . ., 
with parity ±: 

A,± A,±^^.. ^nf -SN^ -A + 1 -A + 3 A-1 
^j'a +2zr/to + 0(e ), a= , 

^ ^ ^ (2.1.3) 

where Imx'^''^ = (mod Zt) and Imx^'~ = t/2 (mod Zt). Complex numbers 
xf'"^ is called a center of a string. Due to Lemma 1.3.2 we may assume that 
Im.x'^''^ = and Ima;^'~ = t/2. We denote the number of ^-strings with parity ± 
byti(A±). 

Note that assumption i) of Theorem 1.4.1 is satisfied for ^-strings with parity 

+, if A = 2£ (mod 2) and for any strings with non-zero real abscissa. Assumption 
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2.2. Ground state and free energy. The result in this section was announced 
in [Take2]. The ground state configuration is specified as follows: u — 0, 
tt(2£,+) = N/2, tt(^,±) = ti(2^,-) = for 2e and centers of 2£-strings dis- 
tribute symmetrically around 0. (Hence ^fl^ x^''^ = 0). This is consistent with 
the result of the XXX type model by Takhtajan [Takh], Babujian [Bab], that of 
the XXZ type by Sogo [So], Kirillov and Reshetikhin [KiR] and that of the XYZ 
model by Baxter [Baxl]. 

Multiplying the Bethe equations (2.1.1) for Xj — a — —£ + 1/2, ... 1/2 

(cf. (2.1.3)), and taking the logarithm, we obtain 

^-1/2 

N ^{xj;2i{a + £)r]t) = 

a=-£+l/2 

N/2 pe-i 2e-i 

k=l \m=l m=0 (2f.2.l) 

where we omit the index (2£, +) of xY'~^ and function $ is defined by: 

^/ N 1, 9ii(x + iLit:it) ^ 
^x; ifit) = -r log -ii) + TT. 2.2.2 

We take the branch of ^> so that $(±1/2; z/x) = ^tt, $(0;z//) = 0. Half integers 
in (2.2.1) specify the branches of logarithm. Applying Takhtajan-Faddeev's 
argument [TF2] to our case, they satisfy — Q^-i — We also assume that 
Xj are ordered by j: Xj > Xj-i. (Note that ^{x;ifj,t) is a decreasing function by 
Lemma C.2.) 

We assume that centers of 2£ strings fill in the interval (—1/2, 1/2) with density 
p{x) in the limit N — > oo, 

p{x), Xj — > X, A?" — > oo. 



N{xj+i - Xj) 



Subtracting (2.2.1) for j from that for j + 1 and taking the limit, we obtain an 
integral equation 

^-1/2 

X ^'{x;2i{a + £)rit) ^ -27vp{x) + 

a=-e+l/2 

nl/2 /2^-l 2^-1 \ 

+ / V $'(x - y; 2imr]t) + V ^'{x - y; 2i{m + l)r]t) p{y) dy. (2.2.3) 

•^-1/2 Vm=l m=0 / 

This equation is easily solved by Fourier expansion. Using formula (C.2), we have 
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Let us compute the eigenvalue for this Bethe vector. The expression of the 
eigenvalue (2.1.2) consists of two terms: 

M 



fj^ Oii{x-Xj;it) 



(2.2.5) 



f^^ Oii{x-Xj;it) 



and the eigenvalue is t{x) = Ai{x) + A2{x). Both Ai and A2 contribute equally to 
t{x) in the thermodynamic limit N ^ 00 for £ > 1/2, since 

1 , Ai , eii(x + 2ier]t;it) 

— loEC — = loe: — h 

N ^As ^ dii(x - 2ier]t : it) 

1 ^ f 9ii{x-Xj -i{2£ + l)r]t-it) 0ii{x - xj - i{2i - l)r]t]it) 
^ TV ^ V 011 (x - Xj + i{2e+ l)r]t; it) ^ O^ix - Xj + i{2e - l)r]t;it) 

N^oo-, 9ii(x + 2i£r]t;it) 

> l0£C h 

^ Onix - 2i£r]t : it) 

r eii{x^y-i{2i+l)r]t;it) O^jx - y - 1(21 - l)r]t;tt) \ ^ 
J-i/2\0ii{x-y + i{2e+l)r]t;it) 9 ii{x - y + i{2£ - l)7jt; it) J ^ 

X p{y) dy 

= - 2TTi. 

In the case of the eight vertex model {£ = 1/2) for A > 0, Ai is dominant in 
magnitude, because 

1 ^ Ai ^ 6ii{x + ir]t;it) 1 9ii{x — Xj — 2ir]t;it) 
iV °^ A^ ~ eii{x -irjt : it) ^ ^ ^ Onix - xj + 2ir)t;it) 

Sttz . isin27rnx (2.2.6) 

> TTIX ^ 



-E 



^ ncosh27rn?7t 

the real part of which is positive. (Recall that x — itX, A > 0.) This is a subtle 
difference between £ = 1/2 and higher spin cases, but the final result does not differ 
much. Namely, in the thermodynamic limit, 

logt(a;) log Ai(a;) 
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Substituting (2.2.4) into this, we obtain the free energy (1.1.5): 

- /3/(A) = (const.) + log ^11 (A + 2er]; r) - 27rt(A - r]){l - Mrj) 

sinh7rnt(l — Mr]) sinh27rnt(A — r]) 

— / : — -. z—T . (2.2.8) 

n smh Trnt cosh ^Trnvt 

n=l 

Here (const.) is an unessential term which does not depend on A. 

2.3. Low-lying excitations and S matrices. As is seen in §2.2, the ground state 
consists of N/ 2 2£-strings filling the Dirac sea. In this section we perturb this Dirac 
sea, slightly changing the string configuration. Since we are interested in the two 
particle states, we choose such configurations that reduce to two particle states of 
models associated to trigonometric and rational R matrices [TF2], [Takh], [KiR], 
[Bab], [So]. 

Let us consider the following configurations: 

(I) +) = N/2 - 2, «(2£ - 1, +) = 1, tt(2£ + 1, +) = 1; 
(II) «(2£, +) = N/2 - 1, ^{2£ - 1, +) = 1, «(1, -) = 1. 

We call the Bethe vectors specified by these data excited state I and II respectively. 
In the higher spin XXX case [Takh] , for example, two particle states are specified 
by similar configurations; one (singlet) is the same as I above, the other (triplet) is 
defined by U(2£, +) = N/2 — 1, tt(2£ — 1,+) = 1. Since one-string with parity — goes 
away to infinity when t tends to oo, we can expect that excited state II reduces 
to the triplet state in the rational limit. As a matter of course, when £ = 1/2, 
{2i — l)-string is absent. Hence the following argument needs to be modified, but 
one obtains results for £ = 1/2 by simply putting £ = 1/2 in formulae for general £. 
We do not mention this modification. 

Excited state I. Now we consider the excited state I. We omit the plus sign desig- 
nating the parity, since all string have parity -|-. Multiplying the Bethe equations 
(2.1.1) for a 2£-string xj = x^^^, a = —£+ 1/2, ...,£ — 1/2 with the center xj^, and 
taking the logarithm, we obtain 

e-i/2 

N ^{xf;2i{a + £)r]t)^27:Qf + 8TT£r]{p + 2Ei) + 

a=-e+l/2 

N/2-2 /2i-l 2e-l 

+ E E *(^f - ^'fc'; 2zmr/t) + J2 ^i^f - 4'; 2^(m + l)r/t) 

k=l \m=l m=0 

2^-3/2 

+ J2 -xl''^;2im7jt)+^xf ~xl'-^;2i{m + l)7jt)) ^'^'^■'^^ 

m=l/2 
2^-1/2 

+ Yl iHxf-xl^+^;2imrjt) + ^xf-x^^+^;2i{m+l)rjt)), 

m=l/2 

where x"^^^ are the centers of the {2£ ± l)-strings. 
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The argument of [TF2] apphed to (2.3.1) imphes that there are N/2 vacancies 
for QY^s. Thus there remain two vacancies (holes) left unoccupied by centers of 
2£-strings. 

We renumber the centers of strings as follows. 

(i) 2£-strings: xj, j = 1, . . . ,N/2, j ^ ji,j2, where Q'j^ and Q'j^ correspond to 
holes. Following the argument in [TF2] again, we assume that Xj > Xji if 

3 > f- 

(ii) 2£ — 1-string: X-. 

(iii) 2£ + 1-string: x^. 

In the thermodynamic limit centers of 2£-strings fill the interval (—1/2,1/2) 
continuously with density pi{x) and two holes at xi = limxj^ and X2 = lima;-,^. 
(We abuse indices.) Subtracting (2.3.1) for j from that for j + 1, we obtain 



£-1/2 

E 

a=-£+l/2 



^ $'(x; 2i{a + i)r]t) = -2tt (pi{x) + -^(5(a; - xi) + 5{x - X2)) \ + 

-"-1-1/2 ^ ^ 

1/2 /2^-l 2t-l \ 

+ / V $'(a; - y; 2imr]t) + V ^'{x - y; 2i{m + l)r]t) pi{y) dy 

•^-1/2 Vm=l m=0 J 

^ 2e-3/2 

+ — ^ {^'{x - x-;2imr)t) + ^'{x - x-;2i{m + l)rit)) (2.3.3) 

m=l/2 
^ 2^-1/2 

+ — ^ {^'{x-x+;2imr)t) + ^'{x-x+;2i{m+l)r)t)), 

m=l/2 

for large N. The solution of this integral equation for pi{x) is 

pi{x) = p{x) + -^{(^ix - xi) + a{x - X2) -\-ijJ-{x - X-) +a;+(a; - a;+)), 

(2.3.4) 

where p{x) is defined above, a{x) and U3±{x) are solutions of the following integral 
equations: Integral equation for a{x): 

2tt(t{x) = -2nd{x) 

1-1/2 /2^-i 2^-1 \ 

+ / ^'{x - y; 2hm]t) +^^'{x-y; 2i{m + l)r]t) a{y) dy, (2.3.5) 

•^-1/2 Vm=l m=0 / 

Integral equation for lo-{x): 
27ru>-{x) = 

.1/2 /2^-i 2^-1 \ 

J f ^ ^'{x - y; 2imr}t) + ^ $'(a; - y; 2i{m + l)r}t) \ 0J-{y) dy 
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Integral equation for u!+{x): 

-1/2 /2^-l 2£-l 



I [ ^ $'(a; - y] 2imr]i) + ^ - y; 2i{m + l)r]i) j 

•^-1/2 \m=l m=0 / 

2^-1/2 

+ ^ ($'(a;;2zmr7t) + $'(a;;2z(m+l)r/t)). (2.3.7) 

m=l/2 

They are easily solved by the Fourier expansion explicitly: 



, , 1 sinh Tmt sinh 27rnr?t 



cos 27mx, 



M sinh7rnt(l — Airj) smh.A'Kni'qt cosh.2TTnr]t 3 g~) 

21-1 ^2smh2nn(2i- l)r]t ^ „ 
a;_(a;) = > — ^cos27rna;, (2.3.9) 

n=l ' 

, ^ 2sinh7r?it(l - 2(2£+ 1)?7) ^ /^..^x 

a;+(a;) = -l-> ^— — ^— — ^ cos 27rna;. (2.3.10) 

^ ^ sinh Tmt (1 -4£?7) ^ ^ 

n=l 

Product of the Bethe equations (2.1.1) for the {2£ — l)-string a;_ + 2ia'qt, a = 
—£ + 1, . . . ,i — 1, gives the equation: 

e+i 

N J2 ^{x-;2i{a + e)r]t) = 27tQ^J-^ -{2£-l)47:r]{i^ + 2Ei) 

N/2 2£-3/2 

+ ^ ^2 i^{x- - Xk;2imr]t) + ^{x- - Xk;2i{m+ l)r]t)) 

k=l,kj^ji,j2m=l/2 (2.3.11) 

2e-i 

+ ($(a;_ — x+; 2imrit) + $(a;_ — x+; 2i(m + 1)^7^)) 

This time there exists only one vacancy for Q^~^ which determines the branch. 
We set Q^~^ = 0. In the thermodynamic limit equation (2.3.11) gives an integral 
equation: 

£+1 . 

$(x_; 2i{a + t)r]t) = -^(2-^ - l)47r?7(zy + 2Si) 

a=-£+l 

1/2 2£-3/2 

+ / V ($(a;_ - y; 2imr]t) + $(a;_ - y; 2i{m + l)r]t))pi{y) dy 
J -1/2 (0 



/2m=l/2 (2.3.12) 
^ 2^-1 

+ — ($(a;_ — x^] 2imr]t) + $(a;_ — Xj^; 2i{m + 1)^7^)). 



m=l 

This equation reduces to 

C)p /'X--X1 , 
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On the other hand, product of the Bethe equations (2.1.1) for the {2£ + l)-string 
x+ + 2im]t, a = —£, gives the equation: 



N J2 *(a^+;2z(Q! + £)r/t) = 27rQ5f^^ - (2£ + l)47rr/(zy + 2Si) 
a=-e+i 

N/2 2£-l/2 

+ ^ ^ {^{x+-Xk;2imr]t) + ^{x+-Xk;2i{m + l)r]t)) 

fc=l,fe^ji,j2m=l/2 (2.3.14) 
2£-l 

+ ^2 (^(^+ ~ x-;2imr]t) + ^{x+ - a;_; 2i{m + l)r]t)) 

m=l 



Again only one vacancy for Q"^^^ which determines the branch exists. We set 
q2^+i _ Q rpj^g integral equation in the thermodynamic limit given by (2.3.14) is 



£ 1 

J2 2^(« + i)vt) = --{2e+ l)A7rr]iu + 2Ei) 



N 

a=-£+l 

1/2 2^-1/2 

+ / ^ {^{x+ - y;2im'r]t) + ^{x+ - y;2i{m + l)r]t))pi{y) dy 
•^-1/2^=1/2 (2.3.15) 

^ 2^-1 

+ — ^ {^{x+ - x-;2imr]t) + ^{x+ - a;_; 2i{m + l)r]t)). 

m=l 



and hence 



f ^ ' a;+(y)dy + a;+-^i^ = (l-2(2£+l)r7)Si- (2^ + 1)^.77 

^ (2.3.16) 



by (2.3.4), (2.2.4), (2.3.8), (2.3.9), (2.3.10) and (C.l). 

Let us denote the solution of the Bethe equations for the ground state by a;^. 
Polarization of the Dirac sea of 2£-strings for excited state I is defined by 
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where x — liniN^ooXj. (Sec [K], [JKM].) Subtracting (2.3.1) from (2.2.1) and 
using the integral equation (2.2.3), one can derive the integral equation for J{x): 

.1/2 /2^-l 2^-1 \ 

-2tcJ{x)+ / ^'{x - y; 2imr]t) +^^'{x-y; 2i{m + l)r]t) J{y) dy 

•^-1/2 \m=l m=0 / 

= - 87T£r]{iy + 2Si) 

2£-3/2 

+ ^ {^{x - x-;2imr]t) + ^{x - x-;2i{m+ l)r]t)) 

m=l/2 

2£-l/2 (2.3.18) 

+ {^{x — x^;2imr]t) + ^{x — x^;2i{m + l)r]t)) 

m=l/2 

(2^-1 2^-1 \ 

^ $'(a; - Xa, 2imr]t) + ^ ^'{x - Xa, 2i{m + l)r]t) J . 
m=l m=0 / 

Thus the polarization is determined as 

J(a:) = ^ J^e^™, (2.3.19) 
nez 



A 2£ - 1 4£ - 1 , 

Jo = ?7(i^ + 2Si) - ^+ + (^1 + ^2), 

_ sinh27rn(2£- l)?7t / _2^i^^_ g-^™! + e-^™^ 

<-'n — : : ^ : ~ I ^ 



(2.3.20) 



27rm sinh 4:7rn£r]t \ 2 cosh 27rnr]t 



^ sinh7rn(l-2(2£+l)r/)t / _ 62^^^+62^^(2.3.21) 

27rm sinh 7rn(l — 4^77) t \ 2cosh27rnr7t /' 

On the other hand, by the definition of the polarization (2.3.17), 
.1/2 

2e / J{x) dx = Ei- {2i + l)x+ - (21 - l)x_ + 2£{xi + X2). 

•^-1/2 (2.3.22) 

Combining (2.3.20) and (2.3.22), we obtain 

x+ - ^i±^ = (1 _ 4£r])Ei - 2£iyr]. (2.3.23) 

Now we determine x± in terms of Xi, X2 regarded as free parameters. Wc have 
derived three equations connecting xi, X2 and x±: (2.3.13), (2.3.16) and (2.3.23). 
From (2.3.13) and (2.3.23) follows 

I \-{y)dy = Q. (2.3.24) 

J —X- +X2 

Thus X- = {xi + X2)/2, since uj-{y) < because of (2.3.9) and Lemma C.2. 

171 /o O 1 £:\ ] o oo^ 
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or, equivalently, 



= J; - (2< + 1)/) , Ei = (2<+1)^, x+ = ^i±^ + 



-a;++a;2 

Equation (2.3.26) shows that is uniquely determined for each ly because of 
Lemma C.2. Constraints (2.3.23), (2.3.26) and (1-4.1) on x+ and Ei are simul- 
taneously satisfied if we put 

kr' xi + X2 , kr' 

y 

where k is an arbitrary integer. Recall that r is assumed to be even. Apparently 
there are infinitely many solutions, but in fact only two Bethe vectors are indepen- 
dent in this series: 

Lemma 2.3.1. For two integers k,k' such that k = k' (mod 2), corresponding 
Bethe vectors are linearly dependent. 

Proof. Let us denote i/ and x+ corresponding to k and k' by {u{k),x+{k)) and 
{i'{k'),x+{k')) respectively. Then 

u{k) - u{k') = -{2£ + l){k- k'Y + ^^^r 

= -{2e+l){k-k'y (modr), 
k — k' 

x+{k)-x+{k') = -^—r'- 

Shifting x+ by one means shifting 2£ + 1 of solutions of Bethe equations by one. 
Lemma is proved by Lemma 1.3.2. □ 

If we take /c = 0, then = 0, Si = and x+ — {xi + X2)/2. We call the Bethe 
vector corresponding to this configuration excited state Iq. 

Since r and r' are coprimc, there exists an integer k such that kr' = 1 (mod r). 
Shifting x+ by an integer, we may assume that x+ — (a;i+a;2)/2 + l/2 with a 
suitable v. (See Lemma 1.3.2.) We call this Bethe vector excited state Ii. 

It seems that there are no other solutions for x+ and u, since the integrality 
conditions (i^ e Z and Theorem 1.4.1) are very strong. 

Excited state II. Now we consider the excited state II. As in the case of the excited 
state I, multiplying the Bethe equations (2.1.1) for a 2£-string Xi = , ol = 

—£ + 1/2, — 1/2 with the center x"^'^ and taking the logarithm, we obtain 

^-1/2 

N $(a;f'+;2z(a + %t) = 27rQf + 87r£r/(zy + 2S„) 

a=-^+l/2 

Ar/2-1 /2^-l 2^-1 

+ *(^f - 2imy]t) + Y $(a;f '+ - '+; 2i{m + \)y]t) 

k=l \m=l m=0 
2^-3/2 (2.3.27) 
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where x^~^''^ is the center of the {2£ — l)-string, {xq + it/2} (xq e M) is the 
one- string with parity — , 

N/2-1 

Sii= (sum of aU x/^) - ^ = 2£ ^ a;f '+ + (2£ - + xj'", 

j=i (2.3.28) 

and function ^{x;iiJ,t) is defined by 

^/ X 1, 9 mix + i tit- it) „ 
* x; = ^ log /i) 2.3.29 

« Ooi[x — int;it) 

By the same argument as for the excited state I, there are N/2 + 1 vacancies for 
Q^^'s. Thus there are again two holes of centers of 2£-strings. 
We renumber the centers of strings as follows. 

(i) 2£-strings: Xj, j = 1, . . . , N/2 + 1^ j ^ ji,j2, where Qjf and Qjf correspond 
to holes. Following the argument in [TF2] again, we assume that Xj ^ Xj' 
if j > j'. The string with its center at xn/2+i will be placed at the zone 
boundary a; = 1/2 in the thermodynamic limit. 

(ii) (2£ — l)-string: X-. 

(iii) 1-string with parity — : xo+ ^. 

As in the previous case, we obtain an integral equation for the density of centers 
of 2£-strings, pii{x), on the interval (—1/2, 1/2) 



^-1/2 

E 



J2 *'(a^f ; 2^(« + i)vt) = -27r (puix) + ^{5ix - xi) + Six - X2))) 
-i+1/2 ^ ^ 

.1/2 /2^-l 2^-1 \ 

+ / *'(^ ~ 2imrit) + V $'(x - y; 2i{m + l)rit) pii{y) dy 

^ 2^-3/2 (2.3.30) 
+ — ^ {^'{x - x-;2imr]t) + ^'(x - x-;2i(m+l)r]t)) 

m=l/2 

+ -^*'(a; - xo; 2i{2£ + l)r]t) + *'(a; - xq; 2i{2£ - l)r]t) 
for large A^. Its solution is 

Puix) = p{x) + ^{crix - xi) + a{x - X2) +u>-{x - X-) + uJo{x - xq)), 

(2.3.31) 

where p{x) ((2.2.3), (2.2.4)), a{x) ((2.3.5), (2.3.8)), U-{x) ((2.3.6), (2.3.9)) are as 
defined before, and u!o{x) is a solution of the following integral equation: 

^1/2 /2i?-l 2£-l \ 

2-KUJo{x) = ^'(x - y; 2imr]t) + ^ ^'{x - y; 2i{m + l)r]t) j ujQ^y) dy 

"^-1/2 \m=l m=0 / 

+^'{x- i{2i + l)rit) + V(x; i{2e - l)rit). (2.3.32) 
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The Bethe equations (2.1.1) for the {2£ — l)-string gives the equation: 

£+1 

N/2+1 2i-3/2 

+ Y Yl ~ ^'^^ 2^™^*) + ^{x- - Xk; 2i{m + 1)77^^)3.34) 

fe=l,7^iij2 m=l/2 

+*(a;_ - xq; 2Ur]t) + *(a;_ - xq; 2i{e - l)r)t)). 

Of 1 

We set Qr"' = 0, since there is only one vacancy. The corresponding integral 
equation in the thermodynamic limit is: 

e+i . 

J2 Hx-; 2z(« + i)vt) = --{21- l)47rr/(i/ + 2E„) 

1/2 2^-3/2 

+ y Y {^{x- -y;'2imr]t) + ^{x- -y;2i{m + l)r]t))pu{y)^^^^^^ 

m=l/2 

1 

+ — (*(a;_ - xq; 2iiT]t) + ^{x- - xq\ 2i{e - l)r]t)). 
This equation reduces to 

-X-+X2 



2£- 1 



/ ' u;_{y)dy + xo - ^^-^ = {1 - 4£rj)Eu - 2£iyrj. (2.3.36) 

J-X-+XO ^ 



as before. 

The Bethe equations (2.1.1) for the 1-string with parity — gives the equation: 
iV*(a;o; 2%ir]t) =2'nQ\~ - 47Tr]{iy + 2Sii) 

N/2 

+ Y {^{xo-Xk;i{2e+l)r]t) + ^{xo-Xk;i{2e-l)r]t)) 

fc=i,#jij2 (2.3.37) 
+*(a;o - x-;2i£r]t) + -^{xq - X-, 2i{£ - l)r]t) 

Again we choose the branch Ql~ = 0. The integral equation is 

*(xo; 2ier]t) = - ^47rr]{iy + 2Eii) 
.1/2 

+ / i^{xo - y; i{2l + l)r]t) +^{xo- y; i{2l - l)r]t))pu{y) dy 
•^-1/2 (2.3.38) 

+ -^(*(a;o - x-;2ier]t) + *(a;o - X-;2i{e - l)rit)). 
This is rewritten as 

fXo—Xl 



f 



uMdv = -2n(2J:u + (2.3.39) 
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for J{x): 

^1/2 /2^-l 21-1 



-27tJ{x)+ I [ ^ - 2/; 2zm?7t) + ^ - y; 2i{m + l)?7t) j J(|/) 

"'-1/2 \m=l m=0 / 



= - ^>^T^r}{v + 2Sii) 

2^-3/2 

+ ^ ($(a; - a;_; 2zm?7t) + $(a; - a;_; 2z(m+ 1)77*)) 

m=l/2 

+^{x - xo\ i{2e + l)r]t) + "^{x - xo; i{2£ - l)r]t) (2.3.40) 



2e-i 21-1 



+ ^ - \;2imr]t) + ^ - \; 2i{m + l)?7t) 

m=l m=0 

/2^-l 2^-1 

- \^2,(^'{x-Xa; 2imr}t) + ^ $'(a; - x^; 2i(m + 1)?/^) j . 



a=l,2 \m=l m=0 

Hence the polarization is 

J(a;) = ^ J^e^'^^"^, (2.3.41) 
X , X 1 2£-l 

(2.3.42) 

_ sinh27rn(2£- 1)77* 



X e 



sinh ATrnirjt 

-2Tvinx 



2 cosh 27Tnr]t sinh 7rnt(l — Mrj) ^ ^ ^^-^ 

^ sinh7rn(l-2(2£+l)?7)t ^ ^ g-2™i ^ g-2™2^ 

2 cosh 27mrit sinh 7m (1 — 4£?7)t 

On the other hand, by the definition of the polarization (2.3.17), 

.1/2 

2e / J{x) dx = Sii - {2£ - 1) X- — xq — £ + 2£{xi + X2). 

•^-1/2 (2.3.44) 

It follows from (2.3.44) and (2.3.42) that 

Xo - ^i±^ = (1 - Mr])En - 2£vr]. (2.3.45) 

From (2.3.36) and (2.3.45) follows the same equation as (2.3.24) and thus X- = 
{xi + X2)/2 as in the case of the excited state I. Equations (2.3.39) and (2.3.45) 
imply 

Lo{y) + dy = -i^r^^. (2.3.46) 

-X0+X2 \ 



\-Ah]J " 1 - Ut] 
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Excited 


2£-strings 


n + l-string 


2£ — l-string 


l-string 


States 


parity + 


parity + 


parity + 


parity — 


lo 


density pi 


x+ = 


X- = 






holes xi, X2 


{xx + X2)I2 


{xi+X2)/2 




II 


density pi 


x^ = 


X- = 






holes xi, X2 


{Xi +X2 + l)/2 


{xi+X2)/2 




IIo 


density pn 




X- = 


Xo = 




holes xi, X2 




{xi+X2)/2 


(Xi +X2)/2 


III 


density pn 




X- = 


Xo = 




holes x\, X2 




{xi+X2)/2 


(xi + a;2 + l)/2 



S matrix. Above we found four excited states with two free parameters xi, X2'- lo, 
II, IIo, III- In the rational limit, t — > oo, rj ^ 0, rjt fixed, the string configuration 
of Iq becomes that of the singlet state of the corresponding spin chain, whereas the 
configurations Ii, IIq, Hi seem to approach to that of the triplet states, since the 
one-string with parity — and the string with abscissa x+ = (xi + X2 + l)/2 goes 
beyond the sight. (Recall that real abscissas of strings are rescaled so that they 
fill the whole real line in the limit.) Hence one might expect that these four Bethe 
vectors give four dimensional space of two physical particle states (spin waves) of 
the corresponding spin chain. In fact for the eight vertex model 

log T (x) I excited state — log T{x)\ ground state 

N 1 . X 2.3.47 

= log T(a; — irjt — xi) + log T(a; — irjt — X2), 

where (excited state) means any one of the excited states lo, Ii, Ho, Hi and 

T^i ■ . v-^ sin27rna; , „ 

logT(x) := -— -mx-iy — , (2.3.48) 

^ ^ ^ 2 ^ ncosh27rnr7t' ^ ^ 

n=l 

(see [JKM] for details of calculations). This means that all conserved quantities 
such as momentum P{x) or energy over the ground states are split into two terms: 

, . sin27rna; 

P{X) = —TTX — } — , 

ncosh27rn?7t 

n=l ' 

and thus we can regard these excited states as two particle states of the XYZ spin 
chain. 

For higher spin cases we have not yet computed fused transfer matrix which 
corresponds to the spin chain with local interaction. From the result of the rational 
and trigonomeric models [Takh] , [KiR] , we conjecture that the momentum and the 
energy of physical particles do not depend on the spin £. Based on this conjecture, 
we calculate the S matrix of two physical particles. 

As is discussed in [DMN] , the S matrix of physical particles (spin waves) could 
depend on the way of calculation for the case of higher spin. In order to make our 
standpoint clear, let us recall the calculation of eigenvalues of the 5" matrix in more 

r\ tit Q 1 1 G fr^l 1 r^Tiri n rr 
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function accumlated should be an integer multiple of 27ri. The main contribution 
comes from the momentum P of the particle as iPN, the free phase: 



iPN = -iNttx -iN^ 



°° sin 27Tnx 



N 



^ ncosh27rnr7t 

n=l 

Note that the right hand side of the above equation is eventually expressed as the 
logarithm of (the right hand side) / (the left hand side) of the Bethe equation for 
the ground state (2.2.1): 

£-1/2 

-iN ^ <^{x\2i{oL^l)j]t) 

N/2 /2£-l 21-1 

+ i ^ ^ ^{x - Xk] 2imr]t) + ^ $(a; - x^; 2i{m + l)?;^) 

k=l \m=l m=0 

£-1/2 

iN <^{x;2i{a + £)r]t) (2.3.49) 

a=-£+l/2 

.1/2 /2£-l 2£-l \ 

+ iN V $(a; - y; 2imr]t) + V $(a; - y; 2i(m + l)r]t) p{y) dy 

-^-1/2 Vm=l m=0 / 

= iPN. 

Hence the ground state can be interpreted as a Dirac sea of non-interacting particles, 
since the momenta of particles are integer multiples of 27r/N because of (2.2.1) and 
the above equation. 

For the excited state, however, the phase shift comes not only from this free phase 
but also from the interaction between physical particles. Because of the periodic 
boundary condition which fixes the total phase shift to an integer multiple of 27ri, 
this means that the calculation of scattering phase shift of a physical particle is 
equivalent to the calculation of 0{1/N) shift of the momentum. In other words, 
the S matrix of physical particles can be calculated by splitting the total phase 
shift, an integer multiple of 27rz, into the free phase iPN of order 0{N) and the 
scattering phase of order 0(1). 

We consider the excited state I first. The total phase shift for the physical 
particle with rapidity xi can be read off from (2.3.1) as follows: 

1-1/2 

-N ^ $(xi;2i(a + £)?7t) + 87r£?7(zy + 2Si) + 

a=-£+l/2 

7V/2-2 /2£-l 2£-l 

+ I] ( I] - ^fe^! 2zmr/t) + ^ ^(xi - xf; 2i{m + l)rjt) 

k=l \m=l m=0 

2£-3/2 

+ J2 mxi-x^J-^;2imr]t) + ^xi-x^_^-^;2i{m + l)r]t)) (2-3.50) 

m=l/2 
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which is equal to an integer multiple of 2ni because of (2.3.1). Subtracting the free 
phase contribution iPN (2.3.49) from the total phase (2.3.50), and taking the limit 
A?" — > oo, we obtain the remainder of order 0(1) and thus we can interpret it as 
the scattering phase shift from the above argument. An explicit expression for the 
eigenvalue of the S matrix for the excited state I is as follows: 

i log(±,Si(a;i - x^)) = -%TT£r]{v + 2Ei) 

,1/2 /2i-l 21-1 



+N I I ^ $(a;i - y] 2imr]t) + ^ $(xi - y; 2i{m + l)r]t) 

•^-1/2 \m=l m=0 / 



X 



X {pi{y) - p{y))dy 

2^-3/2 

+ ^ {^{xi-x_;2imr]t) + ^{xi-x_;2i{m+l)r]t)) (2.3.51) 

m=l/2 

2^-1/2 

+ {(^{xi — Xj^]2irnrit) + (^{xi — Xj^]2i{rn + Vjrjt)), 

m=l/2 

The term —Snir]!/ can be interpreted as an effect from the background or boundary, 
while the rest of the right hand side comes from interaction of pseudo-particles. The 
ambiguity of sign comes from normalizations of asymptotic states. The right hand 
side is computed by integrating (2.3.5), (2.3.6), (2.3.7). The result is: 

i\og{±Sj{x)) = 

E°° / sinh7rnt(l - 4£r/ - 2r/) sinh 7Tnt( Mri - 211) \ 
— r-j 7 7-7-^ T— 1 — — — sin27rna; 
n-i \nsmhTTnt[l — Mrj) cosh27rnrjt nsmhATrnirjt cosh27rnr]t J 

2 sinh 'n:nt{Air) — 2r]) . 

sinh47rn£?7t (2.3.52) 

n=l ' 



00 



Esinh7r?it(2r? — (1 — 4£?7)) . , , 
sinh7rnt(l - 4£?7) ^ ^ 

where £ is or 1 for the excited state Iq or Ii, respectively. The first term of the 
right hand side come from holes {a{x — X2) in pi{x) — p{x) of (2.3.51)), the second 
term from the (2£ — l)-string {u)-{x — X-)) and the last term from the (2£-|-l)-string 
{uj+{x-x+)). 

Computation for the excited state II is the same. The result is 
i\og{±Sn{x)) = 

E°° / sinh7rnt(l — 4£?7 — 2?7) sinh7rnt(4£?7 — 2?7) \ 
I _l_ 1 gin 2tttix 

^-i V'^sinh7rnt(l — 4£r]) cosh.27rnrit nsmhATrni-qt cosh27rn?7t / 

2 sinh nnt{4£r] — 2r]) 



n=l 



— : — - — sm vi/tx /f) Q croN 

sinh47rn£?7t {z.o.oo) 

sinh 2Tmrjt 



sr-^ Slllli ^ /( III i . , , 

-|-7r + Tra; + > — ; -— smTrnia; — £), 
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term from the (2£ — l)-string and the last term from the one-string with parity — 
{uo{x - Xo)). 

We fix the signs left undetermined so that the above S matrix is the permutation 
matrix in the non-interacting limit x — and the excited states Iq reduces to a 
singlet while other three states form a triplet as in the rational and trigonometric 
cases. Then, 

On {^-itri;it{l-Airi)] 



Six)\ 

excited state In 
'S'(a^) I excited state Ii = (x) 
'S'( a;) I excited state IIq = So{x) 
'S'( a;) I excited state IIi — So{x) 



en 


(f +it?7;zt(l 


- Airi)) 


010 


(1 - itr];it{l 


- Air])) 




+itr];it{l 


- 4ir])) 


doi 


(f - itr); it{l 


~ Air])) 


Ooi 


{^ + itr];it{l 


- 4ir])) 




{^-it7]:it{l 


- Airi)) 


^00 


{^ + itr];it{l 


- Mr])) 



(2.3.54) 



where 



So{x) 



-2'jTix 



On -itr];4£itr]) 



Oil {I ^ itr]; Uitri) 
and function S(a;; ]£) is defined by 



^{x;l-Ah])S{x]Ah]). 



(2.3.55) 



§(zAt; ^) = exp f y ""^^^'^"^ 
\n~^ nsmhTrntyU cosh27r: 



sin 27rniXt 



(qV;p^,9^)oo 



ntr] 



i\t\ 



■,p^,q^)oo 



2 477 ' 



1 + ^ 

^ 477 



X 



(2.3.56) 



X 



r.^ + ^) ' (1 + ^) r,. (^) 



n 



where p = e q = e ^'^^^ = p"^. (See Appendix C for definitions of notations. 
The last equality is due to (C.7).) This S factor was found by Freund and Zabrodin 
[FZ]. 

Comparing (2.3.54) with (A. 12), we come to the following conclusion (cf. [FIJKMY]): 

S{x) oc i?(A; it(l - Air])). (2.3.57) 

3. Comments and discussions 

. In this paper we have studied the eigenvectors of transfer matrix of higher spin 
generalizations of the eight vertex model by means of Bethe Ansatz. Apparently 
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two-particle states which would degenerate to a singlet and a triplet in the limit, 
?7 — ^ 0, t — s> cxD. Therefore we can expect that Bethe Ansatz for our case gives as 
many eigenvectors as that for the rational and trigonometric cases. 

• Developments of the theory of quantum afl&ne algebras in the last decade provided 
algebraic tools such as vertex operators and crystal basis for investigation of the 
models associated to trigonometric R matrices. This kind of algebraic method is 
still hard to apply to the models associated to elliptic R matrices because of the 
lack of knowledge on "elliptic affine algebras" which should be an aflfinization of 
the Sklyanin algebra in an appropriate sense. Foda, fohara, Jimbo, Kedem, Miwa 
and Yan [FIJKMY] proposed a candidate of this algebra. In their formulation a 
relation of the type RLL — LLR* plays an important role, where R* is essentially 
the S marix of two particle states of the XYZ model. Their argument was based 
on Smirnov's conjecture, which is supported by the result of the present paper. 

The algebra which Foda et al. propose is considered to be symmetry of the XYZ 
spin chain in a thermodynamic limit. We can also expect that their algebra is also 
a symmetry algebra for higher spin models which we considered in this paper. On 
the other hand, it is still unknown whether finite size models could have symmetry 
of the Sklyanin algebra, since a reasonable coproduct for the Sklyanin algebra is 
not yet found. 

. We considered only such excitations with finite number of holes and finite number 
of strings which have length A, A 2£, even in the thermodynamic limit. This 
is because we wanted to determine the two particle S matrix. But when we want 
to calculate thermodynamical quantities like entropy or specific heat of the model, 
string configurations with non-zero hole density and non-zero densities of A-strings 
(A 7^ 2£) are essential. (See [YY], [TS], [KiR].) Extensive thermodynamics of 
the XXX, XXZ models and their generalization to higher spin cases has quite 
interesting features [R2] , and are also related to dilogarithm identities [Bab] , [KiR] . 
Further study of the thermodynamic Bethe Ansatz for higher spin generalizations 
of the eight vertex model could give deformation of the above features in rational 
and trigonometric models. 
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Appendix A. Review of the Sklyanin algebra 

In this appendix we recall several facts on the Sklyanin algebra and its represen- 
tations from fSkll and \Sk2\ . We use notations in [Ml for theta functions: 
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where r is a complex number such that Im(T) > 0. We denote t — i/r. The Pauh 
matrices are defined as usual: 



1 0\ ,1 _ l^ ,2 _ -A ,3 _ A 

(A.2) 



= 'O li' ll oi' U i' vo -1 



The Sklyanin algebra, Ur,r){sl{2)) is generated by four generators S^, S^, S'^, S^, 
satisfying the following relations: 

Ri2iX - //)Loi(A)Lo2(/x) = Lo2(/x)Loi(A)i?i2(A - /x). (A.3) 

Here A, /x are complex parameters, the L operator, L{X), is defined by 

3 

L(A) = 5^W^i'(A)5«®a«, 

0=0 



6*00(^7; t)' 6*01 (ryiT)' 

i?(A) = i?(A;zt) is Baxter's R matrzx defined by 

3 

R{X) = J2w^Wa-^a-, W^{X) :=W^{X + v). (A.5) 



a=0 



and indices {0,1,2} denote the spaces on which operators act non-trivially: for 
example, 

3 3 

i?12(A) = ® ® ^01 (A) = ^aWS" ® CT« ^ 1. 

0=0 0=0 

The above relation (A.3) contains A and /x as parameters, but the commutation 
relations among S'" (a = 0, . . . , 3) do not depend on them: 

= -iJ^AS^S^U, [S'^,S^].=i[S^,S^]+, (A.6) 

where (a,/?, 7) stands for any cyclic permutation of (1,2,3), [^, = AB ± BA, 
and J«,^ = {W^ - Wl)/iW^ - Wi), i.e., 

'J12 — 



J23 — 

J31 = 



The spin i representation of the Sklyanin algebra, : Ur,r]{sl{2)) — > Endc(0oo) 
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It is easy to see that dim0QQ''~ — 2£ + 1. The generators of the algebra act on this 
space as difference operators: 

(P (S )f)iz) = ^^^^ , (A.8) 

where 

so{z) = 6'ii(?7;t)6'ii(22;;t), 51(2;) = 6'io(?7; r)6'io(22;; r), 

S2{z) = ieoo{r];r)eooi'2z-T), Ssiz) = 6'oi(?7; t)6'oi(22;; r). 

These representations reduce to the usual spin £ representations of U{sl{2)) for 
Ja)3 — >^ (?7 — > 0). In particular, in the case £ — 1/2, S"' are expressed by the Pauli 
matrices a"- as follows: Take (^00(22;; 2t) — ^10(22;; 2t),6oo{2z; 2t) + 6io{2z; 2t)) as 
a basis of Qqq- With respect to this basis have matrix forms 

A/2fca^ Jooim r)Ooi{v; r)eio{r]- T)^ii(r7; r) 
' = '^oo(0;r)^oi(0;r)^,o(0;r) " " ^^"^^ 

Since the relations (A. 6) are homogeneous, an overall constant factor in a represen- 
tation is not essential. 

There are involutivc automorphisms of the Sklyanin algebra Ur,n{sl(2)) found 
by Sklyanin [Sk2]: for a = 1, 2, 3, 

Xa : (5°, S'', S\ S^) ^ (-5°, S"", -S\ -S''), (A.IO) 

where {a,b,c) is a cyclic permutation of (1,2,3). Combining these operators with 
p^, we obtain another representation p^- o Xa of Ur,r]{sl{2)), but there is a unitary 
operator Ua intertwining and p^ o Xa [Sk2] : 

U, :G^S+ 3 f{z) ^{UJ){z) = e-^/ + ^ , 

Us :e^^+ 3 f{z) ^{Usf){z) = e-^e-^(4^+-)/ + ^) , 

and U2 = U3U1. Direct calculations show that Xa{p^{S^)) = U-^p\S^)Ua. Oper- 
ators Ua satisfy the relations: Ul = (-1)^^ UaUi, = {-if^UhUa = Uc- 
Baxter's R matrix (A. 5) is a 4 x 4 matrix proportional to that in [TFl]: 



R{X; it) 



fa{X) d{X)\ 

c(A) 6(A) 

6(A) c(A) 

\d{X) a{X)J 



(A.ll) 



where functions a, b, c, d are defined by 

a(A) = Ci^oi(2ztr/; 2it)9oi{itX; 2it)9ii{itX + 2ztr/; 2it), 
6(A) = Cieii{2itr]; 2it)eoi{itX; 2it)eoi{itX + 2itr]; 2it), 
c(A) = Cieoi{2itr]; 2it)en{itX; 2it)eoi{itX + 2itr]; 2it), 
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Obviously eigenvalues of this matrix are 

Oooi^r -itv;it) , , , , 9qi (H^ - itr];it) 

a A + d(X) = C2 ) -f, ^ , a A - d(X) = C2 \ ^ , 

6(A) + c(A)-C,^^^^^^.^^_.^^, 6(A)-c(A)-C,^-^^p^^, 



where 

^2 



C2 ^26-^*^(^+2") X 

9oo{itr]; it)9oi{itr]; it)9io{itr]; it)9ii{itr]; it) 



Appendix B. Proof of the sum rule 

We prove here the sum rule of Aj's, Theorem 1.4.1. See §1.3 for notations. 
Let us introduce a determinant t^{X) of a (r — 1) x (r — 1) matrix, elements of 
which are defined by: 

(1) {j,j + l)-elements = h{X + 2{j - £)r]); 

(2) (j,i)-elements = t{X + 2jr))- 

(3) {3,3 - l)-elements = /i(A + 2(j + 

(4) other elements are 0, 

where t(A) is the eigenvalue of the transfer matrix T(A) on the Bethe vector 
\l/;^(Ai, . . . , Am) and h{z) = {29ii{z))^ . (This determinant is related to a fused 
model.) Since t(A) is an entire function of A (recall that the transfer matrix T(A) 
itself is an entire function of A), f {X) is an entire function of A. (In other words, 
the analyticity is a consequence of the Bethe equations as noted in §1.3.) Our third 
assumption (see Theorem 1.4.1) is 

iii) t''(A) is not identically zero. 

Let us define t^{X) by 

r(A) := Q{X + 2ri)Q{X + 4r/) . . . Q(A + 2(r - l)r])e{X). (B.l) 

Then because of (1.3.12) 

function t'"(A) is expressed as a determinant of a matrix such that 

(1) (j, i + l)-element = a+(A + 2jri); 

(2) (j,i)-element = a_(A + 2jr]) + a+(A + 237]); 

(3) (j, J — l)-element = a_(A + 2j?7), 

(4) other elements are 0. 
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This determinant can be easily expanded, the result being 

r 

= J] a_(A + 277) . . . a_ (A + 20' - I)r7)a+(A + 2jrj) ...a+{X + 2{r- l)rj) 

i=i 

=h{X + 2{e + l)rj) ...h{X + 2{r-e- l)v)Q{X){fo{X) + ■ ■ ■ + fr-i (A)/,^-^) 

where 

21 fc-1 r 

fk{X)^l[h{X + 2{k-e + j)fj)l[Q{X + 2jfj) n Q(A + 2jr7). 

j=l j=l j=k+2 

By the definition (1.3.11), Q(A) has automorphic property: 
Q(A + l) = (-l)^^-'^g(A), 

g(A + r) = A^^(l+r+2A)-.^r.+2v.^E,^l A,Q^;^^)_ ^ " ^ 

Hence /fc(A + 277) = A+i(A) (/,(A) = /o(A)) and F(A) = /o(A) + • • • + ^_i(A) has 
a period 2r/: F{X + 2r]) = F{X). 
Now wc proceed in four steps. 

Stepl. First we show that Q(A)F(A)/Q(A + 277) . . .Q(A + 2(r - 1)77) is an entire 
function of A. Since 

^ ^ g(A + 277)... g(A + 2(r- 1)77) 

= MA + 2(1 + 1),) ...hiX^2ir-i- lW Q(,^,,)^'.'^ff^',(,_/)y ) 

is an entire function of A, we have only to show that any zero of the denominator is 
notazeroof /i(A + 2(£ + l)77).../i(A + 2(r-£-l)77). Zeros of /i( A) is mod Z + Zr. 
Hence the last statement is true if assumption i) of Theorem 1.4.1 is fulfilled. 

Step2. We show that F{X)/Q{X + 277) . . . Q{X + 2(r - 1)77) is an entire function of 
A. 

As a consequence of Stepl, we know that only possible poles of F(A)/Q(A + 
277) . . .Q{X + 2{r — 1)77) exist at zeros of Q{X). Suppose Xj is a pole of F{X)/Q{X + 
277)... Q(A + 2(r- 1)77). Then 

ordA, F{X) < OTdx, (Q(A + 277) . . . Q(A + 2(r - 1)77)) ^ ordA, F{X) + ordx, Q(A). 

(B.5) 

Here ordA^ is the order of zero at A^-. Assumption ii) of Theorem 1.4.1 says that 
there is an integer a such that ordAj+2ar7 Q(A) = 0. On the other hand periodicity 
F{X + 277) = F{X) implies that ordA,+2a77 F{X) = ordx^ F{X). Therefore 

ordA,+2a^(Q(A + 277) . . . Q(A + 2(r - 1)77)) 

= ordA, (g(A + 2(a + 1)77) . . . Q(A + 2(r + a - 1)77)) 

= ordA, (Q(A + 2a77) . . . Q(A + 2(r + a - 1)77)) 

= ordA, (Q(A + 277) . . . Q(A + 2(r - 1)77)) ^^'^^ 
> ordxj F{X) = oidxj Fix) 



32 



TAKASHI TAKEBE 



Step3. Now we show that even F{X)/{Q{X)Q{X + 2r]) . . . Q{X + 2(r - l)r])) is an 
entire function of A. It follows from Step 2 that for any j — 0,1, . . . ,r — 1 

Q{X + 2jr^)F{X) 
Q{X)Q{X + 2?7) . . . Q{X + 2(r - l)r^) 

is an entire function. Suppose F[X) / {Q{X)Q{X + 2rf) . . . Q(A + 2(r — 1)^7)) has a pole 
at Aq. Then Aq should be a zero of Q{X + ^jrf), j = 0, . . . , r — 1. Taking (B.3) into 
account, this contradicts assumption ii). 

Step4. We have shown that F{X)/G{X) is an entire function where G{X) = 
Q(A)g(A + 2?7) . . . Q(A + 2(r - 1)77). Using (B.3) and 

/i(A + l) = (-l)^/i(A), 

h{X + r) = e-"^^(^+^)-2^^^/t(A), 

we obtain 

^(A + 1) _ F{X) 

G{X + l) G{xy ^ ' ' 

Holomorphy of F/G and periodicity (B.7) make it possible to expand F/G into a 
Fourier series: 



(F/G) (A) = Y.^F/GW-- 



IninX 

Substituting X + t into this expansion and comparing with (B.8), we find that each 
coefficient should satisfy 

{F/G)n = (F/G)ne2"^^("-")^-2^?=i 

Since Imr > 0, there exists only one n — n2 such that {F/G)n2 7^ and it satisfies 
(z/ — n2)T — 2 Yl!jL\ =• ~'"'0 G ^- Putting ni = u — 77,2, we have 2 X^^i A^ = 
no + niT. 

It follows from the above argument that t'^(A) has the following form with a 
suitable integer n: 

e{X) = const. e2'^^^^/i(A + 2{i + l)r]) . . . h{X + 2{r - £ - l)r])Q{Xf. 

(B.9) 

Appendix C. Table of useful functions 
Here we collect properties of functions used in Chapter II. 

Logarithm of quotient of theta functions. A function $ defined by (2.2.2), 

^x; ifit) = log -ii) + TT, 

I dii{x — I lit; it) 
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Hence 

d ^, , ^ sinhTrnfl — 2u)t 
— ^(x;iiJ,t) = -27r { 1 + 2 y ^cos27rna; . 

\ n=l 

A function * defined by (2.3.29), 



^ , . , 1 , 6*01(0; + i\xt\ it) 
* x; tl^t) = -r log -^i) 

has the following Fourier expansion if < < 1/2: 

siiili 2Tvnat 

W(a;;i//t) = 2 > sin27rnx. 

nsinhTrnt 

n=l 

Hence 

d ^ , , ^ sinh 21771111 

— W(a;; = 47r > — — cos 2TTnx. 

dx ^-^ smh Tznt 

n=l 

Lemma C.2. For < a < 6, a series 

y- sinh7rna ^27rmx 
^-^ sinhTrnfe 

zs positive for x gM.. Here the term n = is understood as a/b. 
Proof. Define a function f{y; x) by 

' sinh iryb ' 
f{0;x) = a/b. The Fourier transformation of this function is 

/oo 
f{y;x)e-^-'y^dy 
-00 

2sin(a7r/6) ^-2^\x-i\/b 



> 



6 (e-2,r|x-^|/6 + cos(a7r/6)) + sin2(a7r/6) 
By Poisson's summation formula we have 

E^^— E/(-)-E/(-)>o. 

This proves the lemma. 

For example d/dx{^{x;iiJ,t)) < 0, because of Lemma C.2. 
q-T function. A g-analogue of the F function is defined by 

where {x;q)oc = n^o(l ~ ^^'^)- 

Double infinite product {x;qi,q2)oo is defined by (.x; 52)00 = Tl^i.Tia 
xa?^a^^). lix + v — z + w. the foUowine: relation holds: 
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